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KARMEN, the KArlsruhc Rutherford Medium Energy Neutrino experiment, is lo- 
cated at the pulsed spallation neutron source ISIS of the Rutherford Appleton Lab- 
oratory. In the ISIS beam stop f^, v,, and are produced from the Tr^-zi"*" decay 
chain at rest with energies up to 52.8 MeV. Besides a very low contamination, ISIS 
stands out for its unique time structure. This allows for a highly sensitive search for 
i^/^^i^e-oscillations with the KARMEN detector, a 56t segmented liquid scintillation 
calorimeter with very good time, energy and position resolution. In 1996 an additional 
third veto counter was installed within the 7000 1 steel blockhouse that shields the 
detector against cosmic induced background. Covering the detector from all sides it 
strongly reduces the background of cosmic induced high energy neutrons by a factor 
of 40. Here we present the data taken after this major upgrade from Feb. 1997 until 
Feb. 1999. Since there is no indication for any beam excess, an upper limit for the 
oscillation is deduced using the Unified Approach based on a maximum likeli- 
hood analysis. The result, as all the data presented before, questions the interpretation 
of the LSND beam excess as an indication for i>^j,^Ve oscillation. 
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1 Introduction 

The search for neutrino oscillations is one of the most fascinating topics of modern 
particle physics. The KArlsruhe Rutherford Medium Energy Neutrino experiment 
KARMEN searches for neutrino oscillations in different appearance ( z/^^ z/g and 
i^^J.^ ) and disappearance modes {u^^ ^x^)- The physics program of KARMEN 
also includes the investigation of i/-nucleus interactions as well as the search for 
lepton number violating decays of pions and muons and a test of the V-A structure 
of the /i"*" decay 0. In the following we present the result of the search for v^^Ve 
oscillation on the basis of the data taken from February 1997 until February 1999 
(KARMEN 2 data) after the experiment upgrade in 1996. The data taken before 
the upgrade from 1990 - 1995 (KARMEN 1 data) is not included in the analysis. 
Such a combined analysis would yield a much lower sensitivity due to the relatively 
high cosmic induced neutron background of the KARMEN 1 data. In the data set 
presented here we measure the expected number of background events. Therefore 
we used the Unified Approach based on a maximum likehhood analysis to derive 
a 90 % confidence interval. 

2 Neutrino Production and Detection 

The KARMEN experiment utilizes the neutrinos produced by the neutron spalla- 
tion source ISIS of the Rutherford Appleton Laboratory in Chilton, Oxon, UK. An 
intense beam (200 yuA) of protons is accelerated to an energy of 800 MeV by a rapid 
cycling synchrotron. The two parabolic proton pulses of 100 ns base width and a 
gap of 225 ns are produced with a repetition frequency of 50 Hz (duty cycle is 10~^). 
The protons are stopped in the compact tantalum beam stop. Apart from spallation 
neutrons a large number of pions is produced and stopped immediatly within the 
target. While almost all vr" undergo nuclear capture, the 7r+ decay at rest (DAR) 
into /i"*" and v^. The /i"*" are also stopped within the target and decay at rest via 
/i"*" — > e^ + z/g + v^. The minor fraction of vr" that decays in flight (0.65 % relative 
to 7r+ DAR) with an again suppressed subsequent /i~ decay leads to an extreme 
small z/g contamination of z^e/^'e < 6.2 ■ 10~^ 0. The energy spectra of the neutri- 
nos are well defined due to the DAR of both the 7r+ and The z/^ from 7r+-decay 
is monoenergetic with E(z/^)=29.8 MeV; the continuous energy distributions up to 
52.8 MeV of the z/g and z/^ can be calculated using the V-A theory and show the 
typical Michel shape. Therefore ISIS is a unique, isotropic source of z/^, z/g and z/^ 
from 7r+-/i+ DAR that stands out for its time structure, the small Ve contamination 
and the well defined time and energy distribution of the produced neutrinos. 

These neutrinos are detected with the KARMEN detector, a segmented calorime- 
ter of 56 1 of liquid scintillator. The matrix structure consists of 512 (32 rows x 
IQ columns) optically independent modules with a cross section of 17.4 x 17.8 cm^ 
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and a length of 353 cm. The segmentation is made of thin double acrylic walls sep- 
arated by a small air gap. Every module is read out by two 3 inch photo tubes at 
each end. The position of an event within one module is given by the time dif- 
ference between the photo tubes at both ends. The optimized optical properties 
of the organic liquid scintillator and an active volume of 96% result in an energy 
resolution of as = 11.5%/ ^E[MeV]. Gd203 coated paper within the module walls 
provides an efficient detection of thermal neutrons owing to the very high capture 
cross section of the Gd(n, 7) reaction (cr ^ 49000 barn). The KARMEN electron- 
ics is synchronized to the ISIS proton pulses to an accuracy of 2 ns to fully exploit 
the time structure of the neutrinos. The detector is well protected against beam 
correlated background as well as the hadronic component of the cosmic radiation 
by a blockhouse made of 7000 1 of steel. Cosmic muons entering or stopping close 
to the detector are identified by the two inner veto counters. The innermost veto 
covers the calorimeter from four sides and consists of modules identical to those of 
the calorimeter but half their width. The second veto counter is made of 136 plastic 
scintillator modules that shield the detector from five sides. With this configura- 
tion (KARMEN 1), the dominant background for the search for v^^Ve oscillations 
were high energetic neutrons produced by cosmic muons within the steel blockhouse. 
To eliminate this background source an additional third veto counter made of 136 
plastic szintillator modules with a total area of 300 m^ was installed in 1996 [§. It 
was placed right inside the steel blockhouse such that every muon is detected that 
could produce a neutron within the blockhouse at a distance of up to 1 m from 
the detector. With this configuration (KARMEN 2) cosmic induced background is 
considerably reduced by a factor of 40. 

3 v^^—^Ve oscillation signature 

The probability for v^-^Ve oscillations can be written in a simplified 2 fiavor de- 
scription as 

Am^T 

P(z/^ ve) = sin^ (2 e) ■ sin2(1.27— ^) (1) 

El, 

where L is given in meters, is the neutrino energy in MeV, and Am^ denotes the 
difference of the squared mass eigenvalues Am^ = Iml—m^l in eV^/c"^. The signature 
for the detection of a. Ue is a spatially correlated, delayed coincidence of a positron 
from p ( i^e , e+ ) n with energies up to Ee+ = Ejj^-Q = (52.8-1.8) MeV = 51.0 MeV 
followed by the 7 emission of either of the two neutron capture processes p ( n, 7 ) or 
Gd ( n, 7 ). The p ( n, 7 ) reaction leads to one 7 with an energy of -^(7) = 2.2 MeV 
whereas the Gd ( n, 7 ) process leads to 3 7 in average with a sum energy of 8 MeV. 
The positrons are expected with a 2.2 fj,s exponential decrease of to the /i"*" decay 
after beam on target. The time difference between the positron and the capture 7 
is given by the thermalization and diffusion of the neutron. 
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Figure 1: Signature of sequences of a positron (prompt event) and the correlated gammas from 
the neutron capture reaction (sequential event) that are expected for the D^—fV^ oscillation in the 
KARMEN detector: a) visible energy of the positron for three different values of Am^ as given 
by a Monte Carlo (MC) simulation; b) time of the positron relative to beam on target; c) visible 
energy of the delayed gammas from the nuclear neutron capture on either the free protons in 
the scintillator or the gadolinium in the segmentation; d) time difference of the neutron capture 

reaction relative to the positron. 

To suppress cosmic induced background, a positron candidate is accepted only 
if there is no activity in the central detector and in both inner veto counters up to 
24 /is before. When only the outermost third veto counter was hit, a dead time of 
14 IIS is applied. 

The unique signature of the p ( i/g , e"*" ) n reaction allows already for a strong dis- 
crimination of cosmic and neutrino induced background. The following cuts are 
introduced to maximize the sensitivity of the experiment: The positron has to be 
detected in a time window from 0.6 — 10.6 fis after beam on target with its energy 
in the range from 16 — 50 MeV. The sequential gamma must have an energy below 
8 MeV and has to be correlated in space (within 1.2 m^) and time (5 — 300 /xs) to 




5 



the positron. For these cuts the total detection efficiency is - shghtly depending on 
Am^ - approximately 20%. The expected signature for oscillation sequences in the 
KARMEN detector is shown in Fig. 1. 

Table 1: Expected sequences from background reactions within the cuts specified above. Given 
are the mean values and their errors. Shown in the two last rows are the number of expected 
p (Pe je"*" ) n reactions from the v^—^Vf, oscillation for high Am^ (= 100 eT^^) assuming maximal 
mixing (i.e. sin^ (2 0) = 1) and the number of actually measured sequences. All numbers are 
given for two different energy windows from 16 — 50 MeV and 36 — 50 MeV respectively. 



Background reaction 


events (E>16MeV) 


events (E>36MeV) 


-^^C ( z/e , e~ ) "^^Ng.s. reaction 
V induced random coincidences 
Ve contamination from ISIS 
cosmic induced background 


2.6 ± 0.3 
2.3 ± 0.3 
1.1 ± 0.1 
1.9 ± 0.1 


0.00 ± 0.01 
0.09 ± 0.03 
0.31 ± 0.03 
0.56 ± 0.07 


total expected background 


7.8 ± 0.5 


0.97 ± 0.08 


measured sequences 


8 





p ( i^e ! e"*" ) n reactions for sin^ (2 0)= 1 


1605 ± 176 


712 ± 78 



4 Background Sources 

One of the main advantages for the search of v^^Ve oscillations with the KARMEN 
experiment is that the expected background is not only very small but also known 
with a high percision because most of it can be independently measured by applying 
different cuts. There are only four different sources of background: 

• Pe induced sequences caused by the charged current reaction ^^C ( i/g , e~ ) ^^Ng.s. 
where the subsequent (3 decay of the ^^Ng.s. (r = 15.9 ms) occurs within the 
first 300 yus. 

• Neutrino reactions that have a random coincidence with a low energy event 
from the natural radioactivity inside the detector. 

• The small intrinsic z/g contamination from the 7r^-/i~ decay chain in the ISIS 
target. 

• Undetected cosmic muons which enter the detector or produce high energy 
neutrons via deep inelastic scattering in the inner part of the steel blockhouse. 

The only background source not accessible to direct measurement is the contam- 
ination. It is calculated using a detailed MC simulation of the ISIS target as well 
as all pion and muon production and decay or capture reactions 0. Table 1 lists 
all background reactions and gives the number of expected events as well as their 
errors for the above defined cuts. 
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5 Maximum Likelihood Analysis 

The maximum likelihood (ML) analysis is the most powerful method to infer the 
strength of a possible signal or to derive an upper limit if such a signal is not seen. 
Because of some advantages over other methods we use here the Unified Approach 
[^] recommended by the PDG to derive a 90 % confidence interval from our ML 
analysis. For this ML analysis every background reaction and a possible oscillation 
signal are taken into account with their different probability density functions for 
the time and energy of the prompt event (the e"*") as well as the energy, and the time 
and position difference of the sequential event (the neutron capture) relative to the 
prompt event. The relative contributions of the individual background sources to 
the total number of background sequences is fixed whereas the number of oscillation 
sequences is allowed to vary freely. As an additional information, the likelihood 
function (LF) is weighted with a factor that is the conditional poisson probability 
of the number of inferred background sequences given the expectation value of the 
total background. The resulting LF depends on Am^ and sin^ (2 0) only, and thus 
for a given Am^ only on the number of oscillation sequences No infered (or the 
number of background sequences , for NB=^totai-^o)- 

For the Unified Approach we divided the relevant [Am^;sin^ (2 G)] parameter 
space in the interval [{10~^eV'^ (10~^, 1)] using a logarithmically equidis- 
tant grid of 90 x 72 points [1^. At every point on the grid we generate 8000 MC 
data samples according to the expected background and the given values of Am^ 
and sin^ (2 0) for this point. To these data samples the same ML analysis as to the 
experimental sample is applied. For every MC sample of this specific point on the 
grid one calculates the logarithm of the likelihoodratio A/o(yf[L(Am^; sin^ 20)]mc 
of the value of the LF at its global maximum in the [Am^;sin^ (2©)] parameter 
space to the value of the LF at the given point on the grid for which the sample 
was generated. This procedure gives a characteristic MC generated distribution of 
Alog[L{Am^ ; sin^ 29)]mc for every point on the grid which is then compared to the 
A/o5f[L(Am^; sin^ 20)]exp value of the experimental data set (i.e. the logarithm of 
the ratio of the experimental LF at its global maximum to its value at a given point 
on the grid). The 90% confidence interval (C.I.) is the set of points on the grid 
for which the experimental value A/o(7[L(Am^; sin^ 2G)]exp is smaller than at least 
10% of aU MC generated A/o^[L(Am2; sin^ 2G)]mc (see Fig. 2). If, on the other 
hand, for given parameters [Am^;sin^ (2 0)] A/o(y'[L(Am^; sin^ 2G)]exp hes in the 
upper 10 % tail of the MC distribution this point on the grid does not belong to the 
90% confidence interval. The upper 90% confidence limit (C.L.) as shown in Fig. 4 
is the upper limit of the 90% C.L. The interpretation of this C.L. is that for all 
parameter combinations Am^ and sin^ (2 B) on this curve 90 % of a large number of 
hypothetical KARMEN experiments would have seen a larger "signal" (i.e. smaller 
Alog[L{Am^; sin^ 20)]exp) than the one actually observed if - and this is important 
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- the true parameters of the y^^^e oscillation were in this region of the parameter 
space. 
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Figure 2: Monte Carlo generated distribution of the logarithm of the hkelihoodratios 
A;o.g[L(Ani2;sin^2e)]MC for Am^^ lOOeV^ and sin2(2e)= 0.001. The experimental value 
A/of;[L(Am^; sin^ 20)]exp for this point in the parameter space is to the left of the upper 10% 
tail of the MC distribution. Therefore this parameter combination belongs to the 90% C.I.. 



6 Results and Conclusion 



The results presented here are based on the data recorded in the measuring period 
from February 1997 to February 1999 which corresponds to 4670 C protons on target. 
Within the cuts defined in Sect. 3 we find 8 sequences as shown in Fig. 3. Since we 
expect a total background of 7.8 ± 0.5 sequences there is absolutely no indication for 
a beam excess. For this data set the above described analysis leads to a 90 % CL. of 
sin2(2e)=2.1-10-3forlarge Am2 (i.e. 100 eV^). The90%C.L. as a function of Am^ 
can be seen in Fig. 4. Also shown is the sensitivity of the KARMEN experiment. 
The sensitivity of an experiment is defined as the mean confidence limit a large 
number of identical experiments would yield if there was no oscillation. The actual 
limit is slightly "better" then the sensitivity with sin^(2G)=2.3 ■ 10^^ for large 
Am^. Naively one would expect that the sensitivity is at slightly lower sin^ (2 6) 
than the exclusion curve and not vice versa because we measure more events (0.2) 
than the expected background. This apparent contradiction is explained by the fact 
that there are no events above 36 MeV where - depending on Am^ - roughly half 
of all oscillation sequences should be and only 0.97 ± 0.08 background events are 
expected (see Tab.l). This leads also to an LF that has its global maximum in the 
region of small negative sin^ (2 0) which - for a null result - is as possible as a global 
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Figure 3: Distribution of tiie expected baclcground sequences for tlie prompt energy (a) and time 
(b) distribution as well as the sequential energy (c) and time (d) distribution. Also shown are the 
8 measured sequences which agree nicely in their shape with the expected background. 



maximum in the region of positive sin^ (2 0). The 90 % CL. from our ML analysis is 
compared in Fig. 4 to the LSND resuh |jll|. It excludes most of the LSND favoured 
region and is thus - as is the fact that there is no event above 36 MeV - strongly 
questioning the interpretation that the LSND beam excess is an indication for v^-^Ue 
oscillations. Furthermore one has to keep in mind that the KARMEN limit that was 
derived from the Feb. 97 - Apr. 98 data set and that is also shown in Fig. 4 puts an 
even stronger constraint on the LSND favoured region. In this data set (which had 
slightly different cuts) we expected 2.9 ± 0.1 background sequences and measured no 
event at all. This yields, using again the unified approach, a 90 % CL. of 1.3 ■ 10~^ 
for large Am^ and a sensitivity of 5.4 ■ 10~^, respectively |T^, |TB[. Although this 



limit was criticized by some people, it is of course valid and correct: If one accepts 
the ansatz of the Unified Approach one must NOT argue that this ( "lucky" ) limit is 
not trustworthy for this would mean to reduce the Unified Approach to absurdity. 
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Figure 4: KARMEN 2 90 % confidence limits according to the Unified Approach compared to 
other experiments: The full line is the 90%C.L. of the data presented here, the dotted line the 
corresponding sensitivity and the dashed line the 90% C.L. derived from the Feb. 97- Apr. 98 data. 
Also shown are the 90% C.L. of the two reactor experiments Chooz [14 and Bugey |1^] and the 
favoured region for i/^^i^e oscillations as reported by the LSND experiment (TTj. Areas to the right 
of the 90 % C.L. are excluded with a probability of more than 90 %. For the LSND result the "99 % 
favoured region" (total shaded area) and the "90 % favoured region" (light-shaded area) are given. 
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